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Abstract 

Going beyond a recently proposed microscopic model [1] for the incommensurate transition in the spin-Peierls 
TiOX (X=C1, Br) compounds, in the present work we start by studying the thermodynamics of the model with 
XY spins and adiabatic phonons. We find that the system enters in an incommensurate phase by a first order 
transition at a low temperature T c i. At a higher temperature T C 2 a continuous transition to a uniform phase is 
found. Furthermore, we study the magnetism in the incommensurate phase by Density Matrix Renormalization 
Group (DMRG) calculations on a ID Heisenberg model where the exchange is modulated by the incommensurate 
atomic position pattern. When the wave vector q of the modulation is near 7T, we find local magnetized zones (LMZ) 
in which spins abandon their singlets as a result of the domain walls induced by the modulated distortion. When q 
moves far away enough from ir, the LMZ disappear and the system develops incommensurate magnetic correlations 
induced by the structure. We discuss the relevance of this result regarding previous and future experiments in TiOCl. 
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1. Introduction 

The interest on spin-Peierls systems [2] has been re- 
newed with the recent discovery of a dimerization in 
the TiOX (X = CI, Br) compounds family. In them, the 
magnetically active Ti ions connected by O ones are 
arranged in a bilayer structure. Ti ion positions on a 
layer are shifted with respect to the other neighboring 
layer, forming something like an anisotropic triangu- 
lar structure. Initially postulated for a resonating va- 
lence bond state [3] , TiOCl turned out to be mainly a 
one- dimensional magnetic system [4] with the Ti d xy 
orbitals pointing towards each other in the crystallo- 
graphic b direction [5] . The high temperature magnetic 
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susceptibility is well described by the Bonner-Fisher 
curve, indicating a nearest neighbor magnetic exchange 
J ~ 660K [5]. The peculiar ingredient in the TiOCl 
phase diagram is the appearance of an incommensurate 
intermediate phase between the high temperature uni- 
form phase and the low T dimerized phase [5,6] com- 
monly found in a regular spin-Peierls system. The tran- 
sition temperatures were experimentally found to be 
Tei ~ 66K and T c2 ~ 92K. In Ref. [6] a very large en- 
ergy gap of about 430A" in the low temperature phase 
and a pseudo spin gap below 135K have been also re- 
ported. 

The origin of the intermediate phase remains contro- 
versial and not well understood yet. The shifted posi- 
tions of a Ti ion in a chain between the two neighbor- 
ing chains, make it plausible to entail that some type 
of competition between the in-chain and out-of-chain 
interactions could be the origin of the incommensurate 
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phase. 

Some of us recently proposed a microscopic mecha- 
nism for the incommensurate phase [1] , where the an- 
tiferromagnetic Ti chains are immersed in the phonon 
bath of the bilayer structure. By using the Cross-Fisher 
theory [7], it was shown that the geometrically frus- 
trated character of the lattice is responsible for the 
structural instability, leading the chains to an incom- 
mensurate phase without an applied magnetic field. For 
TiOCl, our results showed consistency with the tem- 
perature dependence of the phonon frequencies and the 
value of the incommensuration vector at the transition 
temperature. Moreover, we found that the dynamical 
structure factor shows a progressive softening of an in- 
commensurate phonon near the zone boundary as the 
temperature decreases, along with a broadening of the 
peak. These features are in agreement with the exper- 
imental inelastic X-ray measurements [8] . 

The purpose of this work is double. As a first step 
we present the phase diagram as a function of the tem- 
perature of the model proposed in Ref. [1]. On the 
other hand, we theoretically analyze the evolution of 
the magnetism inside the incommensurate phase by 
solving a one-dimensional Heisenberg model with a 
distortion-dependent exchange, using an incommensu- 
rate wave as the modulation of the Ti positions. This 
study is performed using Density Matrix Renormaliza- 
tion Group (DMRG) calculations. 



2. The model and its phase diagram 

The model proposed in Ref. [1] contains the relevant 
magnetic interaction, the phonons and the spin-phonon 
coupling. It includes the Ti atoms on the bilayer struc- 
ture which interact by harmonic forces. For simplicity, 
we consider ionic movements only in the direction of 
the magnetic chains. As the measured magnetic sus- 
ceptibility of TiOCl is well reproduced by a ID Heisen- 
berg model, the assumption that the Ti atoms are mag- 
netically coupled only along the b direction was taken 
into account. Moreover, as a direct exchange seems to 
be the dominant Ti-Ti interaction, it is reasonable to 
consider that J will be modulated by the movement 
of nearest neighbor Ti ions in the chain direction. Our 
spin-phonon Hamiltonian reads: 

H — H p h + H sp h, (1) 

H P h = + X^^pK.J ~ w *+!j) 2 

. Kinter r, \2 . / \2l -i 

H — [{uij - Uij+i) + (Ui,j - Ui+ij-i) J }, 

Hsph — ^ ^ Ji,j^i,j " Sj+lJ, (2) 

i,3 



where mj is the displacement of the atom i of the chain 
j along the direction b of the magnetic chains and the 
Pi j are their respective momenta. Sjj are spin-i op- 
erators with exchange constant Jij = J[l + a(uij — 
Ui+i,j)] along the b-axis of a non-deformed underly- 
ing lattice. K in and Ki n t er are the harmonic force con- 
stants. 

The Hamiltonian (1) contains variables of different 
dimensionality: whereas the spins are one-dimensional, 
phonons are two-dimensional. An exact solution to this 
model at finite temperature is a formidable task, there- 
fore we make two approximations to obtain the phase 
diagram. We consider only the XY coupling in the 
spin-spin interaction and phonons are treated in the 
adiabatic approximation. Spin variables are converted 
to spinless fermions by means of the Jordan- Wigner 
transformation. The resulting pseudofermion-phonon 
model is solved in a finite system by diagonalizing the 
one-electron problem in the background of a deformed 
lattice (let us call each eigenvalue A 7 , 7 = l,...,N a , 
where N a is the the number of sites in each chain) . 

Now we consider deformations of the form: 

= mo cos(q ■ - tt/4), (3) 

where the wave vector q = (0, q) is constrained to the 
chain direction and Rjj is the position vector in the 
two-dimensional Bravais lattice. We have chosen the 
7r/4 phase to have a continuous energy at q = it. The 
total free energy per site is: 

/=^^-£Eln(l + e--n (4) 
s 7=1 

where T is the temperature in units of J and uj 2 (q) = 
£ [K in sm 2 (% )+if lnter (l-cos(^)cos(f ))] is the 
dispersion relation of the two-dimensional phonons. 

We have minimized / with respect to Mo and q at dif- 
ferent temperatures. The results are shown in Fig. 1. 
At low temperatures a dimerized phase is found. This 
is the usual result for an isolated spin-Peierls chain. At 
the first transition temperature T c \ , uo jumps discon- 
tinuously and q moves away from it. This is a structural 
incommensurate phase which arises as a competition 
between the tendency to the dimerization of each chain 
and the frustration in the elastic interchain coupling. 
The value of q remains constant up to a second temper- 
ature T C 2 where uo vanishes and the system becomes 
uniform. The sequence of the two transitions and the 
first order character of the first transition [8,9] exper- 
imentally found in TiOCl can therefore be accounted 
by our model. However the precise fitting of the param- 
eters is a quite complex task which would need the in- 
clusion of the z-component of the spin interaction and 
the non adiabatic effects of the phonons [10]. 



2 




0.06 0.08 0.1 0.12 0.14 

T 



Fig. 1. Phase diagram in the plane ua,T of the model 1 as 
discussed in the text. Parameters are a — 3 and — ^t^r _ ^ g 
We identify each phase with the value of the wave vector of 
the pattern that minimize the free energy. 

3. Magnetic structure in the incommensurate 
phase 

The feature that makes the compound under study 
different from an usual spin-Peierls system is the ap- 
pearance of an incommensurate phase without an ap- 
plied magnetic field. It is therefore quite important 
to characterize the magnetism in this phase. As no 
magnetic interchain interaction was included in the 
model (1), we will consider a one-dimensional Heisen- 
berg model as given by (2) with deformation fixed by 
(3). The model will be solved by DMRG. We will focus 
on the local magnetization and spin-spin correlation 
functions. 

As a reference state to discuss the results, we take the 
dimerized phase because the structural deformation (3) 
with q near n can be thought as a smooth modulation 
of the dimer pattern. In a dimerized phase the spins 
are paired in singlets and the local magnetization is 
zero. There is a gap in the spectrum and the magnetic 
correlation declines exponentially. When the pattern 
(3) is included with qi = 7r(l — to have an integer 
number of wavelengths in the chain, the dimerization 
phase changes its sign 21 times as we travel along the 
chain from one extreme to the other (see the gray line 
of Fig. 2). When I is small enough, the situation has 
some similarities with the domain walls (or solitons) 
appearing in the presence of a high magnetic field in 
regular spin-Peierls systems [11]. We expect that for 
each zero of (3) free LMZ should appear. This is in fact 
the situation we found in Fig. 2a, b and c where the 
mean value < Sf > is shown. We see the appearance of 
LMZ in all the zeros of the deformation except in the 
first, which is greatly affected by the proximity to the 
edge. These LMZ are seen as clouds of polarized spins 
emerging from the singlet background. Moreover, they 
are uncorrelated among them because their distances 
are longer than the correlation length of the dimerized 
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Fig. 2. DMRG results for a lattice of N 3 = 200. We show in 
black the atomic positions for different qi — 7r(l — -^-), with 
I — 1,2,3,4 from top to bottom, and uq — 0.053. Superim- 
posed, we plot the < > to show the appearance of the LMZ. 

state. They are randomly orientated. We expect that 
the magnetic susceptibility will be affected by these 
LMZ. It should increase compared to the dimerized 
phase. This feature has been experimentally observed 
in TiOCl [8]. 

When qi moves further away from tv, the zeros of 
(3) fall within a distance of the order of the correla- 
tion length and the LMZ disappear as seen in Fig. 2d. 
The effect of the structural incommensuration on the 
magnetism is more subtle and a careful analysis of the 
magnetic correlation function is necessary. 

In Fig. 3 we show the spin-spin correlation function 
C(| i - j |) =|< Si ■ Sj > - < Si >< Sj >| and its 
Fourier transform (the static magnetic structure fac- 
tor S(q)). For the three q where the LMZ appear, the 
correlation decreases exponentially to zero. Moreover, 
it falls a bit faster the farther q is from n. This seems to 
be a finite size effect because the modulation interferes 
with the decrease of the spin correlation of the dimer- 
ized zones when a zero appears. So we expect that in 
the thermodynamic limit the dimerized region will be 
large enough for the correlation to decrease as in the 
dimerized phase. 

A more dramatic change in the behavior of 
C(| i-j\) is found for q<t where the LMZ do not 
appear. In this case the correlation decreases more 
slowly than in the previous cases. For comparison we 
show the result of a uniform Heisenberg model (UHM) 
where the quasi-long range antiferromagnetic order is 
expected. Note that C(| i — j ) for q<t approaches the 
one of the UHM, a quite different behavior than that 
of q\, 52 and 53. We will undertake a careful analysis of 
this long distance behavior for different modulations 
in a forthcoming paper. 

More importantly, there is a modulation of the long 
distance behavior with half the wavelength of that 
which has the imposed structure. The magnetic struc- 
ture can no longer be interpreted in terms of a weak 
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Fig. 3. (a) DMRG results for the spin-spin correlation functions 
C(l * ~~ j I) f° r deformations with the qi as given in Fig. 2. We 
also show the correlation of a UHM for comparison. In the inset 
we show a log-log plot of C{\ i — j ). (b) Fourier transforms 
of the correlation functions of (a). 



magnetic behavior inside the incommensurate phase 
from a structure with LMZ to another without LMZ 
but with an incommensurate correlation function. This 
should be checked in future magnetic experiments. 

In summary we have found that a system with 
one-dimensional antiferromagnetic chains coupled to 
two-dimensional phonons in an anisotropic triangular 
lattice accounts for the appearance of the incom- 
mensurate phase in TiOCl. The magnetism in this 
intermediate state includes uncorrelated LMZ near 
the dimerized phase, and incommensurate correlations 
with no LMZ towards the uniform phase. In the last 
case, the modulation found in the correlation function 
is related to the modulation of the imposed distortion 
by a factor of two in their wavelengths. 

This work was supported by ANPCyT (PICT 1647) , 
Argentina. 



perturbation of the dimerized state. Instead, the mag- 
netism becomes incommensurate as a response to the 
modulated exchange. In some sense, the situation is 
the dual to the one recently studied in some multifer- 
roic materials [12]. There, an helical magnetic struc- 
ture induces a lattice distortion which finally produces 
the ferroelectricity, whereas in our case the lattice dis- 
tortion is the driving force which produces an incom- 
mensuration of the magnetic state. 

These features are even more apparent in the struc- 
ture factor S(q) shown in Fig. 3b. For qi, qi and q$ 
there is a single maximum at q ma x = tt whose intensity 
decreases as q moves away from 7r. This maximum is 
rounded as in a pure dimerized chain [13]. It can be fit- 
ted near q = n by the law ^ q _ 7r y\ + ^-i (corresponding 
to a free massive boson propagator [14]) with £ ~ 3.3. 
This asserts our previous statement that imposing 
distortions with wavector qi near enough to n, the 
long range magnetic correlation behaves as in a dimer- 
ized chain, decreasing exponentially with a correlation 
lenght of the order of £. For qi = qi the behavior 
changes qualitatively. Two additional maxima appear 
at one side and the other of the highest peak at n, cor- 
responding to the long distance magnetic modulation 
found in C(| i - j |). Moreover, the peak at the center 
is not longer rounded but it behaves similarly to the 
UHM. This signals a tendency towards a quasi-long 
range incommensurate magnetic correlation of the 
system. Let us make a brief discussion about the ap- 
plication of these results to TiOCl. As measured from 
X-ray experiments [8,9], the modulation wave vector 
is q ~ 3.05 at T C 2 and it moves slightly to n when de- 
creasing the temperature towards T c \ where it jumps 
discontinuously to 7r. So we expect a change of the 
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